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We propose an effective two-dimensional quantum non-linear sigma model combined with classical 
percolation theory to study the magnetic properties of site diluted layered quantum antiferromagnets 
like I^Cui-^M^CU (M=Zn, Mg). We calculate the staggered magnetization at zero temperature, 
M 3 (x), the magnetic correlation length, £(x,T), the NMR relaxation rate, l/Ti(x,T), and the Neel 
temperature, Tjsr(x), in the renormalized classical regime. Due to quantum fluctuations we find a 
quantum critical point (QCP) at x c ~ 0.305 at lower doping than the two-dimensional percolation 
threshold x v ~ 0.4f . We compare our results with the available experimental data. 
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The discovery of high temperature superconductivity 
in La2- x Sr a; Cu04 has motivated an enormous number of 
experimental and theoretical studies of this and related 
materials. La2Cu04 has attracted a lot of interest be- 
cause it is a classical example of a quantum Heisenberg 
antiferromagnet (QHAF). I^CuCU is a layered quasi- 
two-dimensional (2D) QHAF, with an intraplanar cou- 
pling constant J (J/fcs « 1500 K) much larger than the 
interplanar coupling J± (« 10~ 5 J) [jjj. The quantum 
nonlinear sigma model (QNLctM) is probably the sim- 
plest continuum model with correct symmetry and spin- 
wave spectrum that reproduces the low-energy behavior 
of a QHAF. It has been successfully used [Q to explain 
many magnetic properties of La2- x Sr x Cu04 

In this paper we propose a QNLerM allied to classi- 
cal percolation theory to study the site dilution effect 
in La2Cui_ ;c M; C 04, where M is a non-magnetic atom. 
While the theory of disordered classical magnetic systems 
is fairly developed || we still lack deep understanding of 
the behavior of the site diluted QHAF Q . As we show 
below the interplay between quantum fluctuations and 
disorder leads to new effects which cannot be found in 
classical magnets. In particular we show that long-range 
order (LRO) is lost before the system reaches the clas- 
sical percolation threshold. Furthermore, we have only 
two independent parameters in the theory: the spin- wave 
velocity cq (rs 0.74 eV A/h ||) and the bare coupling 
constant go (~ 0.685 of the clean system (x = 0). 
The results for the staggered magnetization, correlation 
length, NMR relaxation rate and Neel temperature are 
derived without any further adjustable parameters. 

Our starting point is the 2D site diluted nearest- 
neighbor isotropic Heisenberg model 



H= jJ2p(r i )p(r j )S i -S j 



(1) 



where p(r) is the distribution function for Cu sites: 
p(r) = 1 on Cu sites and p(r) = on M sites. Al- 
though translational invariance has been lost in (|l|), the 
Hamiltonian retains the SU(2) invariance for rotations 



in spin space. Since the symmetry is continuous Gold- 
stone's theorem predicts the existence of a gapless mode 
in the broken symmetry phase. The ordered phase is 
characterized by a finite expectation value of the magne- 
tization, n = (S(Q)), at the antiferromagnetic ordering 
vector Q = (n / ao , n / ao) (ao = 3.8 A). 

In the pure system, in accord with the Hohenberg- 
Mcrmim- Wagner theorem, LRO for a system with con- 
tinuous symmetry is only possible at finite temperatures 
in dimensions larger than 2. In the absence of disorder 
the system has a Goldstone mode which is a spin wave 
around Q with energy E(k) and linear dispersion rela- 
tion with the wave- vector k: E(k) = ?ic|k|, where c is 
the spin-wave velocity. This dispersion relation is a con- 
sequence of the Lorentz invariance of the system. In the 
paramagnetic phase, where the continuous symmetry is 
recovered, all excitations are gapped because order is only 
retained in a region of size £. In this case the excitations 
have dispersion 



E(k) = Hc^k 2 + l/£ 2 . 



(2) 



Now consider the case where quenched disorder is 
present. Spin- wave theory, which can only be applied 
to ([!]) at T = 0, predicts that Lorentz invariance is 
lost even for an infinitesimal amount of impurities [Q. 
The dispersion changes to k ln(fc) and the spin waves be- 
come damped at a rate proportional to k when k — > 0. 
These results (strictly valid in 2D and T = 0) are not 
directly applicable to the systems in question which or- 
der at finite temperature g. At finite temperatures and 
weak disorder we can consider the criterion established 
by Harris for the relevance of disorder in critical phe- 
nomena Firstly, we can classify the phase diagram 
of the pure system as [0: renormalized classical (RC) 
where £,(T) diverges as exp(To/T) (where To is a char- 
acteristic temperature scale - see (pt)); quantum criti- 
cal (QC) where £(T) oc 1/T; quantum disordered (QD) 
where £(T) « £o is constant. If we imagine the pure 
system being divided into regions of size £, each part 
will have fluctuations in the microscopic coupling con- 
stant (5, say) which by the central limit theorem are 
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proportional to the square root of the number of spins 
oc £ 2 in that region. That is, there are statisti- 
cal fluctuations of order 6g(£) oc l/y/N(£) oc l/£. On 
the other hand the thermal fluctuations in the system 
are of order 5T(£) oc l/ln(£/a ) in RC, ao/£ in the QC 
and vanishingly small in the QD region. For the critical 
behavior of the system with weak disorder to be essen- 
tially the same as for the pure system one must require 
that 6T(£) > 5g(£) when £ > a . Observe that this 
condition is always fullfilled in the RC regime and there- 
fore we expect the critical behavior to be the same as 
in the pure system, that is, described by a QNLcrM 0. 
In the QC and QD regimes the situation is not clear 
because 5T(£) ~ <5<?(£) and therefore the effect of dis- 
order is strong. We conjecture that in these regimes 
the critical behavior is different from the one described 
by a QNLcrM. In this work we focus entirely in the RC 
regime. Having these results in mind we can apply clas- 
sical percolation theory to ( l|) |H],[n| • The main param- 
eters of the problem depend on geometrical factors such 
as the probability of finding a spin in the infinite clus- 
ter Pqq(x) (ps 1 — x, for x <C 1) and the bond dilution 
factor [jiJI A(x) (~ 1 — irx + nx 2 /2) (in the expressions 
below Poo(x) and A(x) are valid for all x as given by 
the numerical simulations fllJl ). In the classical case the 
spin stiffness p s {x) is related to the undoped stiffness 
by p s {x) — A(x)p s (0), while the transverse susceptibil- 
ity is given by \±(x) = (Poo{x)/A(x))x±{0) so that Q 
p s (x) = c 2 (x)x±(x). 

In this paper we propose an effective field theory which 
is valid for Tjv < T < J/ks and combines the Lorentz 
invariance implied in (Q), the Harris criterion and the 
results of percolation theory. In percolation theory, be- 
sides the infinite cluster, we always have finite clusters. 
A finite cluster of size L has discrete energy levels and 
therefore a gap of order he/ L. In what follows we assume 
£ 3> L and ignore the contribution of finite clusters to the 
magnetic properties and focus entirely on the physics of 
the infinite cluster. It is obvious from the definition of 
p(r) that on average (p(r)) = P o0 (x). Furthermore, site 
dilution implies that n 2 (r) = p(r). Thus, on average we 
have Q (n 2 (r)) = P 00 (a;). In the continuum limit of (|l|) 
the Harris criterion discussed above indicates that in the 
long- wavelength low-energy limit the magnetic properties 
of the site diluted problem can be described in terms of 
an effective QNLcrM: 

Z = J DnS [n 2 - P^x)] exp{-S eff /h} , 



where 



lit, 



S eff = 1/2 / dr I dr X x(x) \d T n\ 2 + p s (x) |Vn| 2 (3) 



and r is the imaginary time direction with [3 — 1/ (fc^T). 
Equation (Q) leads to a natural description of the un- 
doped system and provides an effective field theory for 
the QNLcrM in the presence of impurities. Moreover, it 



has incorporated the correct properties of the classical 
percolation problem added to the quantum fluctuations 
of the QHAF. It is very simple to show by a change of 
variables that the action in (fcf) can be rewritten as 



S, 



eft 
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2g(x) 



dr / cZr(d M n) 



(4) 



where g{x) = hc(x)/p s (x) is the effective coupling con- 
stant of the theory. Moreover, because of the contin- 
uum limit the theory has an intrinsic ultraviolet cut-off 
A(x) = 2^/7rP oc (x)/ao which is fixed by the total num- 
ber of states. In writing (^) we have not included the 
topological term. In a random system one suspects that 
this term vanishes as in the pure 2D case jl5| . Neverthe- 
less there are always statistical fluctuations in a random 
system which are of order y/Nj, where Ni is the number 
of M ions. Thus, the topological term has importance as 
we discuss at the end of the paper. 

The great advantage of (|4|) is its simplicity and close re- 
lationship to the description of the undoped problem. In 
this paper we use the large N approach for the QNLcrM 
which has been so successful in describing the undoped 
system [^9|. At zero temperature, a critical value of 
the coupling constant g c (x) separates the RC from the 
QD region. g c (x) — 4ttP oc (x)/ A(x) can be obtained 
from the saddle-point equation for |l9) . The ratio of 
the coupling constant to the critical coupling constant is 
g(x) = g{x)/g c (x) = g / 'Paa{x), which implies that non- 
magnetic doping drives the system from RC region to QD 
region at x c where Poo(x c ) = g at T = 0. The critical 
concentration x c is completely determined by the value of 
g a in the undoped case. Using the dilute result for Poo (a;) 
and g — 0.685 we find x c ~ 0.3 which is indeed smaller 
than the percolation threshold x p rs 0.41 p0| . This result 
has to be contrasted with classical calculations p(J where 
long range order is lost at percolation threshold only. We 
also performed a one-loop renormalization group analy- 
sis and calculated the zero temp erature st aggered mag- 
netization Q M s (x) = M (x)y/l -g{x). Here M (x) 
is the classical staggered magnetization for the perfect 
Neel spin alignment and the remaining factor is due to 
quantum fluctuations. Thus, the local average magnetic 
moment is given by 



(n(x)} M s (x)/M (x) _ l-g(x) 



(MO)) M s (0)/M Q (0) 



1-5(0) 



(5) 



Observe that the average local moment indeed vanishes 
at x c . For the undoped case, @ predicts that the max- 
imum measured magnetic moment of Cu ion is 0.56/ig 
which agrees with the measured value 0.6 ± 0.15/xb p2| . 
It is also in good agreement with the existing experimen- 
tal sublattice magnetization measured by /xSR for various 
doping concentrations as shown in Fig. [l]. Notice that 
for the Ising magnet p{x) only deviates from p(0) at x p . 
The larger reduction of the moment in the QHAF is due 
to quantum fluctuations present in the QNLcrM. 
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The magnetic correlation length £ can be directly cal- 
culated from the QNLcrM. The int erpola tion formula 
from the RC to the QC region reads 



Z(x,T) = 



ehc(x) \ exp (2tt p R , s (x)/k B T) 
4irp RtS (x) + k B T 



(6) 



where pr :S (x) — p s (x)[l — ~g(x)] is the renormalized spin 
stiffness. This result agrees very well with the Monte 
Carlo simulations in a large temperature range in the 
undoped case |24[| . As far as we know the only existing 
neutron scattering results for magnetic correlation length 
are for the pure system and x = 0.05 [jl6). In Fig. || we 
plot the available data and the prediction of our model 
given in (^|). As it is well known, samples with x = 0.05 
have problems with the Oxygen stoichiometry []l6|| . Ex- 
cess O introduces mobile holes in the plane which produce 
strong frustration effects which are not accounted for in 
our theory. Thus, direct comparison between the theory 
and experiment for this sample is problematic, especially 
at high temperatures. Thus, only new experiments with 
controlled O content can directly test our theory. 

Chakravarty and Orbach j25[ have calculated the nu- 
clear spin-lattice relaxation rate of Cu for La2Cu04 us- 
ing the dynamical structure factor from the QNLcrM. A 
detailed calculation was done in Refs. 18 2q|. These cal- 
culation can be easily extended for the doped case. Here 
we just quote the result for A£ ^> 1: 



1 



Ti(x,T) 



= 1 2 P OQ (x)V2^S(S+l) 



x e(Aj.-4Poo(a;)B) Jl 



2A ± B 
A\ + 4B 2 



[(A± - 4P 0O (z)B) £ 2 + AP^Ba 2 In (£A)] 



3w e (x)£a (ln(CA)) 2 

where 7 is the nuclear gyromagnetic ratio, A} 
and B = 83 kG are the hyperfine constants llj 



= 80 kG 
and 



U e {x) = A(x)\ 



! / 2J 2 k 2 B zS(S 



1) 



3?T 



(where z is the number of nearest neighbor spins) is the 
corrected Heisenberg exchange frequency. Fig. H shows 
the NMR relaxation rate normalized to the high tem- 
perature value as given by the experimental data and 
the result of our calculations. The growth of the relax- 
ation rate at low temperatures is due to fast growth of 
£. As the system approaches the QCP one starts to see 
the crossover from RC to the QC regime where the £ 
grows like 1/T leading to slower growth of the relaxation 
rate. This behavior is clearly seen in the data since for 
x = 0.11 where growth is very slow from 800K down to 
400K. The agreement between data and theory is again 
quite reasonable. 

The 3D Neel order can be obtained from the weak in- 
terplane coupling Jj_ and it is given by H,EUp7J : 



«() 



which is a transcendental equation for Tn(x). The inter- 
planar coupling constant is insensitive p^j to doping be- 
cause the change in lattice parameters is negligible [ p9[ . 
In the undoped case the Neel temperature Tjv(0) is of 
order of 315 K. The initial suppression rate of the Neel 
temperature with doping, / = — d\n(T]sr(x))/dx, when 
x — > can be directly obtained from (R) and, due to 
quantum fluctuations it is much faster than in the Ising 
case (dashed line on Fig.[|). We find I f=s 4.7 in good 
agreement with the data. Indeed, in Fig. ^ we show our 
theoretical results in comparison with various different 
experimental measurements. The critical concentration 
x c for which the system loses long-range order by moving 
from the RC region to the QD region is approximately 
0.305, in agreement with the loss of long-range order at 
zero temperature as given in (|s|). Finally, it is also easy 
to show using the procedure given in ref. |j0| that the 
topological term will lead to induced moments close to 
the impurities. These moments interact through a ran- 
dom magnetic exchange of order Je _ (°° a: V£( a: > T ). This 
effect can lead to order of the induced moments in the 
paramagnetic phase, as seen experimentally [j3l"| , j32|| . 

In conclusion, we have proposed an effective QNLcrM 
to describe the magnetic diluted QHAF. Our model com- 
bines the result of classical percolation theory and the 
quantum fluctuations of the Heisenberg model. Although 
our model is fairly simple it gives a good quantitative 
description of the magnetism in La2Cui_^M ;c 04. The 
success of our model in describing the physics of the RC 
regime is due to the fact that the 2D correlations are very 
long at finite temperatures and the effect of disorder in 
the critical behavior is rather weak. Disorder induces 
quantum fluctuations in the system which lead to the fi- 
nal destruction of LRO at x c . This effect is not found 
in classical magnets where LRO is solely determined by 
the percolation problem. Finally, our arguments indi- 
cate that a new approach is required in the QC and QD 
regions where the NLctM is probably not applicable. 
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x (Zn concentration) 

FIG. 1. Effective moment at T — as function of x (nor- 
malized relative to the undoped case) and the experimental 
data for La2Cui-2,Zna;04 101. 
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FIG. 2. Inverse correlation length as function of tempera- 
ture for x = (solid line) and x = 0.05 (dashed line). The 
open (x = 0) and solid symbols (x = 0.05) are the neutron 
scattering data Jl(| . 
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FIG. 3. The Cu 1/Ti normalized to the undoped case at 
high temperatures (T = 900 K). The lines from top to bottom 
are for x = (solid), x — 0.025 (long dash), x = 0.08 (short 
dash), and x — 0.11 (dotted dash), x — (solid circles - NQR 
data and x — 0.025 (open tri ang le), x — 0.08 (open 

square) x = 0.11 (cross) (NQR data 0). 
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FIG. 4. Neel temperature normalized to the undoped case. 
Solid line: [jj; dashed line: Ising result. Experimental data for 
La2Cui-iZna;04: Solid and open circles are NQR and /iSR 
data, respectively [Jl8| ; straight and up-side-down triangles are 
magnetic susceptibility data crosses are magnetization 
data [||. 
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